Hypergeometric representation of a four-loop vacuum bubble 
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In this article, we present a new analytic result for a certain single-mass-scale four-loop vacuum (bubble) 
integral. We also discuss its systematic e-expansion in d — 4 — 2e as well as d — 3 — 2e dimensions, the coefficients 
of which are expressible in terms of harmonic sums at infinity. 



1. Introduction 

Recently, a number of calculations have 
reached the four- loop level p^. In most cases, such 
high-loop computations proceed by first reducing 
the (typically very large number of) momentum- 
space integrals to a small number of so-called 
master integrals, a step that can be highly au- 
tomatized 2 . In a second step, these master in- 
tegrals have to be computed in an e-expansion 
around the space-time dimension d of interest. 

This second step is much less amenable to au- 
tomatization, at least if one is interested in ana- 
lytic expressions for the coefficients. Indeed, very 
few four-loop master integrals are known analyt- 
ically. It is the purpose of this paper to add one 
master integral to this set. 

As a basic building block, many of the above- 
mentioned calculations use single-mass-scale vac- 
uum bubbles. These arise e.g. when there are 
large scale hierarchies in the physics problem, 
such that asymptotic expansions can be used. 
Having no external momenta, and propagators 
with masses or m only, members of this class of 
integrals are functions of the space-time dimen- 
sion d only (for convenience, we set m — 1 in the 
remainder of this paper), hence the coefficients of 
their Laurent expansions in e are pure numbers. 

Define the (dimcnsionlcss, but d-dependent) 1-, 
3- and 4-loop integrals 
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At x = 1, all three examples are master integrals. 
While J(l) and B2(l) can be computed analyti- 
cally in terms of Gamma functions, we would like 
to know T 3 (l). 

2. Difference equation 

Using a slight generalization of the abovemen- 
tioned reduction step, one can systematically de- 
rive sets of difference equations for the integrals. 
For an integral f{x), they have the generic form 

YliLo c if( x + *) ~ 9( x id), with polynomial co- 
efficients Ci(d), and where N is the order of the 
difference equation and g(x) its inhomogeneity, 
which is supposed to be known analytically. 

First-order difference equations can be solved 
directly in terms of one initial value. 

In the homogeneous case, this can be done triv- 
ially in terms of Pochhammer symbols (a) = 
T (a + x) /Y (a): 
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For these two cases, we know the initial values 
J(l) and B 2 (l) analytically in d dimensions. The 
value of J(l) depends on the choice of integration 
measure, which we do not need to specify here 1 , 
and a simple integration gives 
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r(}-i)r(i- 



The 4-loop integral T3 is an example for the 
inhomogcneous case, 
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This difference equation has been solved numer- 
ically in 3d [3] and 4d 0]. We will solve this 
equation analytically in the following section. 

3. Solution 

One can immediately write down a solution of 
Eq. (0 in terms of a sum, as can be easily seen 
by doing the first couple of iterations, 

T a (x + 1) = T 3 (x ) ( f[ c(i) j + 
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In this case however, it is the initial value T 3 (l) 
that we would like to compute. What we instead 
know about the integral T 3 (x) is its behavior at 
the boundary, T 3 (x 3> 1) oc J(x ^> 1). 
To proceed, let us rewrite Eq. JfjJ as 
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7r d / 2 r(i - d/2). 



T 3 (z+1) 



T{x + l)T(x + 9 - 3d) 
r(a; + 3-d)r(a; + 6-2d) 

r(j + i)r(j + 9 - 3d) 
r(j+3-d)r(i+6-2d) 



r(a;o + 2- d)T(x + 5 - 2d) 



(8) 



T(x )T(x + 8 - 3d) 

The left-hand side does not depend on x, so 
the same has to be true for the right-hand side. 
Hence, we can freely choose the value of x at 
which we want to evaluate the right-hand side; 
in particular we can choose 1 > 1 or even per- 
form the limit lim^^oo. 

Using Stirling's formula T(x + a)/T(x + b) x = 
x a-b a _|_ q f x -i\\ we gee fag^ ^he g rs ^ term on 

the right-hand side of Eq. © behaves like x 1 ^/ 2 
at large x, so it vanishes for d > 2. So we are left 
to evaluate 
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Multiplying by (1) -/(l) ., the sum is just a gen- 
eralized hypergeometric function j+\Fj at unit 
argument, 
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xo + 6 — 2d, «4 = .to + 1 — I and «5 = xo + 3 — d. 
From this, setting xo = 1, we finally get 
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We will now turn to expanding this function 
around d = 4 and d = 3. 

4. Expansion in 4d 

Following the notation of Ref. 0|, we define 
_ T 3 (l) 



T3 " [^(l)] 4 ■ 
For d = 4 - 2e, Eq. (JTUJl then reads 
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An e-expansion of the 3F2 can be achieved via 
Algorithm A of Ref. . We have used the pack- 
age XSummer [7] to expand in e, in terms of har- 
monic polylogarithms H(z) jH] of unit argument. 
In a second step, we have used the package Sum- 
mer to rewrite the H(l) in terms of a minimal 
set of numbers, which are equivalent to harmonic 
sums S'm(n) at infinity |10|. 

We have coded both expansions in FORM 11 , 
and obtain (note the absence of an e 2 term) 
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where £„ = C( n ) are values of the Riemann Zeta 
function, and s$ a — 5*5,3(00) « 1.041785 is the 
only non-zeta value appearing up to this order. 
Up to the e 5 term, these coefficients have already 
been obtained (by fitting high-precision numeri- 
cal values) in Ref. 0|. Higher-order coefficients 
are new. To 30 digit accuracy, Eq. i|13|) reads 

T 3 = + 0.250000000000000000000000000000 

+ 0.500000000000000000000000000000 e 

- 0.186630129462637144901701950176 e 3 

- 116.322454283613146131176255284 e 4 

- 37.9603995633682022478604992948 e 5 

- 3339.84136081794624883976306535 e 6 

- 580.540498497887663112519574990 e 7 

- 68729.9111560055991193326976118 e 8 

- 7274.24628251698563686391222230 e 9 

+ 0(e w ) . 

5. Expansion in 3d 

For d = 3 - 2e, Eq. 0UJ) reads 
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XSummer cannot expand the 3 F 2 around half- 
integer indices, so unfortunately we cannot pro- 
ceed the same way as above. There are algorithms 
for expansion around rationals p/q |12| . which 
work only if rationals are balanced between nu- 
merator and denominator. The case at hand is 
unbalanced. 

To balance it, we may make use of the Euler 
identity for Gaufi' hypergeometric function 

which allows us to rewrite 
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This gives 6 transforms, 4 of which are balanced. 
Choosing one of those possibilities 2 , 
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where (n) i = are Pochhammer symbols. 

The sum is of "type A" according to the classi- 
fication in ^21; a class that can be expanded in 
terms of harmonic sums, as we will show now. 

Using (e), = -YZA-lT + for two 
terms in the numerator, (ae) i and (6e) i; and solv- 
ing one of the two new sums, the above sum reads 
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2 Note that both generalized hypergeometric functions we 
need can be conveniently summarized introducing only one 
parameter, a = 2 and a = 4, respectively. 



The next step is to expand the Pochhammer 
symbols in e, using 12 
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It becomes clear now why the sum had to be bal- 
anced: the Pochhammer symbols (|) ■ and (1)^ 
cancel between numerator and denominator. Ex- 
panding the exponentials and rewriting products 
of S-sums as single S-sums leaves us with sums of 
generic form J2i S%(2i)/i n to be solved. 

To this end, we use the trick employed in |12j . 
introducing a delta function on the integers: 
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Putting it all together, and doing some trivial 
algebra with Gamma functions, 
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For expanding the Gamma functions around in- 
teger and half-integer arguments, we use 
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